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We report on the first measurement of the rapidity distribution da/dy over nearly the entire 
kinematic region of rapidity for e + e~ pairs in the Z-boson region of 66 < Af ee < 116 GeV/c 
and at higher mass M ee > 116 GeV/c 2 . The data sample consists of 108 pb _1 of pp collisions at 
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^J~s = 1.8 TeV taken by the Collider Detector at Fermilab during 1992-1995. The total cross section 
in the Z-boson region is measured to be 252 ±11 pb. The measured total cross section and da/dy 
are compared with quantum chromodynamics calculations in leading and higher orders. 

PACS numbers: 13.85.Qk, 12.38.Qk 
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Most measurements at high energy proton-antiproton 
colliders are performed in the central rapidity produc- 
tion region, \y\ < 1. A model dependent extrapolation 
for \y\ > 1 is needed to extract the total cross section for 
hard processes such as top quark production or W and 
Z boson production. This extrapolation is made using 
Monte Carlo programs (e.g. PYTHIA [[[)), which incor- 
porate quantum chromodynamics (QCD) calculations in 
leading order (LO) or next to leading order (NLO). A pre- 
vious measurement of the rapidity distribution, da/dy, 
for dimuon pairs in the Z-bosons mass region was lim- 
ited to \y\ < 1 In this communication, we present the 
first measurement of da/dy for e + e~ pairs in the Z-boson 
mass and high mass region over nearly the entire kine- 
matic region of rapidity. At the Tevatron pp collider, the 
kinematic limit at the Z-boson mass is \y\ = 3.0, while 
we measure \y\ up to 2.8. The da/dy distributions are 
compared to the predictions of QCD in LO and NLO. 
This measurement is also relevant for precision W boson 
mass measurements at hadron colliders, where TV's are 
reconstructed using ev and \iv pairs from the Droll- Yan 
process. 

In hadron-hadron collisions at high energies, massive 
e + e~ pairs are produced via the Drell-Yan || process. In 
the standard model, quark-antiquark annihilation form 
an intermediate 7* or Z (7* / Z) vector boson, which then 
decays into an e + e~ pair. In LO, the momentum fraction 
X\ (#2) of the partons in the proton (antiproton) are re- 
lated to the rapidity M, y, of the boson via the equation 
X1.2 = (M I \fs)e ±y . Here s is the center of mass energy, 
and M is the mass of the dilepton pair. Therefore, dilep- 
ton pairs which are produced at large rapidity originate 
from events in which one parton is at large x and another 
parton is at very small x. Since the quark distributions 
for x up to 0.9 are well constrained by the deep-inelastic 
lepton scattering experiments j^j, comparisons of data 
and theory for da/dy, and the total cross sections pro- 
vide a test of the theory, e.g. missing NNLO or power 
correction || terms. 

The e + e" pairs are from 108 pb _1 of pp collisions at 
\fs = 1.8 TeV taken by the Collider Detector at Fer- 
milab j§ (CDF) during 1992-1993 (18.7 ± 0.7 pb" 1 ) 
and 1994-1995 (89.1 ± 3.7 pb" 1 ). CDF is a solenoidal 
magnetic spectrometer surrounded by projective-tower- 
geometry calorimeters and outer muon detectors. Only 
detector components used in this measurement are de- 
scribed here. Charged particle momenta and direc- 
tions are measured by the spectrometer, which consists 
of a 1.4 T axial magnetic field, an 84-layer cylindrical 
drift chamber (CTC), an inner vertex tracking cham- 
ber (VTX), and a silicon vertex detector (SVX). The 
polar coverage of the CTC tracking is |r/| < 1.2. The 
pp collision point along the beam line {Z VCItcx ) is deter- 
mined using tracks in the VTX. The energies and direc- 
tions Q| of electrons, photons, and jets are measured by 
three separate calorimeters covering three regions: cen- 
tral (|?7| < 1.1), end plug (1.1 < |?7| < 2.4), and forward 
(2.2 < \rj\ < 4.2). Each region has an electromagnetic 



(EM) and hadronic (HAD) calorimeter. 

In a previous letter, we presented da/dP^ of Z bo- 
son B. This analysis is an extension of the da/dPr mea- 
surement. The da/dPr analysis has three categories of 
e + e~ pairs: central-central (CC), central-end plug (CP), 
and central- forward (CF). This analysis extends the sam- 
ple to the forward rapidity region by including plug-plug 
(PP) and plug-forward (PF) events. The inclusion of 
these events increases the event sample by 20% and al- 
lows for measurement of Z bosons with \y\ up to 2.8. An 
improvement in this analysis is the additional VTX track- 
ing requirements for plug and forward electrons. The 
VTX covers the entire rapidity range in this study, and 
plays an important role in removing background in the 
high rj region which is not covered by the CTC. 

The sample of e + e~ events was collected by a three- 
level online trigger that required an electron in either 
the central or the plug calorimeter. The offline analy- 
sis selects events with two or more electron candidates. 
Since the electrons from the Drell-Yan process are typi- 
cally isolated, both electrons are required to be isolated 
from any other activity in the calorimeters. Electrons in 
the central, end plug, and forward regions are required to 
be within the fiducial area of the calorimeters and have 
a minimum Et of 22, 20, and 15 GeV, respectively. To 
improve the purity of the sample, electron identification 
cuts are applied M. For CC, CP, and CF events, the 
central electron (or one of them if there are two) is re- 
quired to pass strict criteria. The criteria on the other 
electron are looser. A central electron must have a CTC 
track that extrapolates to the electron's shower clusters 
in the EM calorimeter. These clusters must have EM-like 
transverse shower profiles. The track momentum and the 
EM shower energy must be consistent with one another. 
The track is also used to determine the position and di- 
rection of the central electron. The fraction of energy 
in the HAD calorimeter towers behind the EM shower is 
required to be consistent with that expected for an EM 
shower (-Ehad/^em)- The plug electrons must also have 
an EM-like transverse shower profile. The end plug and 
forward electrons are required to pass the Shad / Eem 
requirement and to have a track in the VTX which orig- 
inates from the same vertex as the other electron and 
points to the position of the electromagnetic cluster in 
the calorimeter. The ratio of found to expected hits in 
the VTX is required to be greater than 70% and 50% 
for plug and forward electrons, respectively. The VTX 
tracking efficiency is (97.8 ± 0.3)% for plug electrons and 
(97.0 ±0.9)% for forward electrons. This requirement on 
plug and forward electrons reduce the background rates 
from 11% to 2% for PP events, and from 22% to 4% for 
PF events. 

The data sample is divided into two mass bins: the 
Z region (66 < M ee < 116 GeV/c 2 ) and the high mass 
region (M ee > 116 GeV/c 2 ). After all cuts, the numbers 
of CC, CP, CF, PP, and PF events in the Z mass region 
are 2894, 3811, 621, 1236, and 589, respectively. The 
backgrounds are low and are estimated using the data. 
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The backgrounds in the CP, CF, PP, and PF topologies 
are dominated by jets and are 28 ± 5, 13 ± 3, 24 ± 2, and 
23 ± 3 events, respectively. Because of the CTC track- 
ing requirement, the jet background for the CC sample 
is negligible. The CC background is mainly from e + e~ 
pairs from W + W~, t + t~, cc, bb, and ti sources. This 
background, estimated using e ± /i =F pairs (To), is 3 ± 2 
events. In the high mass bin, the numbers of CC, CP, 
CF, PP, and PF events are 61, 59, 9, 18, and 5, respec- 
tively. The mean mass of the high mass bin is 152.6 
GeV/c 2 . 

The acceptance for Drell-Yan e + e~ pairs is obtained 
using the Monte Carlo event generator, PYTHIA [Q, and 
CDF detector simulation programs. PYTHIA generates 
the LO QCD interaction (q + q — > ^*/Z), simulates 
initial state QCD radiation via its parton shower algo- 
rithms, and generates the decay, ^*/Z — > e + e~. To 
approximate higher order QCD corrections to the LO 
mass distribution, a "X-factor" |ll[] is used as an event 
weight: K{M 2 ) = 1 + |(1 + §7r 2 )a s (Af 2 )/27r, where a s 
is the two loop QCD coupling. This factor improves the 
agreement between the NLO and LO mass spectra. (For 
M > 50 GeV/c 2 , 1.25 < K < 1.36.) The CTEQ3L @ 
PDFs are used in the acceptance calculations. Final state 
QED radiation |b| from the J*/Z — > e + e~ vertex is 
added by the PHOTOS @ Monte Carlo program. Gen- 
erated events are processed by CDF detector simulation 
programs and are reconstructed as data. The calorime- 
try energy scales and resolutions used in the detector 
simulation are extracted from the data, as are the cut ef- 
ficiencies and corresponding errors. Simulated events are 
accepted if after the reconstruction they pass the e + e~ 
pair mass, the detector fiducial, and kinematic (Et) cuts. 

In the da/dy measurement, various samples are com- 
bined and binned in y. The da/dy is calculated with 

da _ AN 

dy C Ay A- t-K ' 

The AN is the background-subtracted event count in 
a y bin, C is a bin centering correction, Ay is the bin 
width, the sum r is over the 1992-1993 and 1994-1995 
runs, C r is the integrated luminosity, and eA r is the run's 
combined event selection efficiency and acceptance. The 
backgrounds subtracted from the event count are pre- 
dicted using the data and background samples. The fac- 
tor C corrects the average value of the cross section in the 
bin to its bin center value. Acceptances are calculated 
separately for CC, CP, CF, PP, and PF pairs. They are 
combined with the corresponding event selection efficien- 
cies to give eA r . Figure [l] shows the eA r for events in 
the Z region as a function of y for the CC sample, the 
CC+CP+CF sample, and the CC+CP+CF+PP+PF 
sample, respectively. It indicates that the PP+PF events 
extend the acceptance of the y measurement to \y\ = 2.8. 
We significantly increase our statistics by extending the 
acceptance beyond y > 1.2. The da/dy of e + e _ events 
in the Z mass region is shown in Table 0. 
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FIG. 1. The efficiency times acceptance for e + e~ pairs in 
the Z boson mass region: (a) for the CC sample (asterisks), 
(b) for the CC+CP+CF sample (open circles), and (c) for the 
CC+CP+CF+PP+PF sample (solid circles). 

The systematic errors considered are from variations 
in the background estimates using different methods, the 
background in the efficiency sample, the uncertainty in 
energy resolution of the calorimeter, the choice of PDFs, 
and the distribution of Z used in the Monte Carlo 
event generator. The systematic error from the calorime- 
ter resolution is 0.2% in the low y region and increases to 
0.7% at \y | > 2.5. The systematic error from variations 
in the Z P T is 0.5% at \y\ = and is 2.0% at \y\ = 2.8. 
The systematic error from the choice of PDFs is less than 
1% in \y\ < 2.0 and increases to 2% at \y\ = 2.8. For the 
total cross section measurement, the combined system- 
atic error is 0.6% (excluding the luminosity uncertainty). 
The pp collision luminosity is derived with CDF's beam- 
beam cross section, ctbbc = 51.15+1.60 mb p^JIq |. The 
luminosity error of 3.9% contains the <tbbc error and un- 
certainties specific to running conditions. 

Figures ||(a) and |^(b) compare the measured da/dy's 
to theoretical predictions in the Z mass and high mass 
regions, respectively. The top horizontal axes on these 
figures are the corresponding values of the x\ and X2 
as a function ofy. The predictions are LO calculations 
with CTEQ5L @ PDFs and NLO @ calculations with 
MRST99 Q and CTEQ5M-1 PDFs. The predic- 
tions in Figure |^(a) have been normalized by a factor "F" , 
the ratio of measured total cross section to the predic- 
tion (F=1.51, 1.14, and 1.13 for the CTEQ5L, MRST99 
and CTEQ5M-1 PDFs, respectively). The predictions 
in Figure |2|(b) are normalized to the data using the fac- 
tor F from the Z mass region. We compare the data to 
the theory using statistical errors only. As the x 2 values 
listed in Figure ||(a) indicate, the LO calculation using 
recent LO PDFs does not fit the shape as well as the 
NLO calculation with the most recent NLO PDFs. 

Model independent measurements of the total produc- 
tion cross sections for e + e _ pairs are extracted by inte- 
grating the measured values of da/dy. Because there 
are no data for \y\ > 2.8, we use a NLO calculation 
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with the CTEQ5M-1 PDFs to estimate j*/Z produc- 
tion in that region. The cross section in the region of 
\y\ > 2.8 is about 0.02% of the predicted total cross 
section. The extracted cross section in the Z mass re- 
gion is 252.1 ± 3.9 (stat.) ± 1.6 (syst.) ± 9.8 (lum.) pb. 
The corresponding cr(pp — > Z) ■ Br(Z — > ee) is 253 ± 
4(stat. + syst.) ± 10(lum.) pb. This measurement is in 
agreement with our previous measurements in the di- 
electron [§ (248 ± 5(stat. + syst.) ± 10(lum.) pb, using 
only the CC+CP+CF e + e~ sample), and dimuon g] 
(237 ± 9(stat. + syst.) ± 9(lum.) pb using a CC sam- 
ple) channels. These previous measurements use a QCD 
model to correct for the missing events at high ra- 
pidity. The combined e + e _ and cross section 
is 250 ± 4(stat.+syst.)±10 pb. Since the pp inelastic 
cross section used by CDF in luminosity calculations dif- 
fers from DQ's by +5.9% H|, measured cross sections 
must be renormalized before comparisons to other ex- 
periments. The D0 measurement [^o| of 221 ± 11 pb 
when renormalized to the CDF luminosity measurement 
is 234 ± 12 pb. 
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FIG. 2. da/dy distribution of e + e~ pairs: (a) in the Z 
boson mass (66 < M ee < 116 GeV/c 2 ) region, (b) in the high 
mass (Mee > 116 GeV/c ) region. The M used to obtain xi 
and X2 in (b) is the mean mass over the bin. The error bars 
on the data include statistical errors only. The theoretical 
predictions have been normalized to the data in the Z boson 
mass region. 



The total cross section measurements can also be com- 
pared to QCD calculations. Fixed order QCD calcula- 
tions have uncertainties from PDF measurements and 
corrections from higher orders of perturbation theory, 
i.e., the if-factor. The NLO-to-LO total cross section 
correction is significant: K ~ 1.4. In contrast, the NLO 
total cross section is lower than NNLO |(| prediction by 
only 2.3%. The NNLO prediction with the latest NLO 
MRST99 PDFs is 227±9 pb, where the 4% error is mostly 
from uncertainties Jll| in the NLO PDFs. Although a full 
set of NNLO PDFs is not available, recent estimates JjJ of 
NNLO PDFs indicate that the NNLO PDFs will increase 
the theoretical cross sections by 5%. Given these uncer- 
tainties, the theoretical expectation is consistent with the 
Z cross section measurements. 

The measurement of the Drell-Yan total cross sec- 
tion in the high mass region is 4.0 ± 0.4 (stat. + syst.) ± 
0.2 (lum.) pb. The corresponding prediction of the total 
cross section from the NNLO QCD theory using MRST99 
PDFs is 3.3 pb. 

In summary, the rapidity distributions of e + e~ pairs in 
the Z boson mass and high mass region have been mea- 
sured for the first time over nearly the entire kinematic 
region. This measurement uses a new tracking technique 
in the high rapidity region to reduce the background and 
uncertainties associated with it. In addition, unlike the 
previous measurement of the total cross section, this mea- 
surement is model independent. 
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TABLE I. da/dy distribution of e + e~ events in the mass 
range 66 < M ee < 116 GeV/c 2 . The first and second errors 
are statistical and systematic, respectively. The 3.9% lumi- 
nosity error is not included. Here y is the bin center value. 



y 


da/dy [pb] 


y 


da/dy [pb] 


0.05 


70.78 ±3.27 ±0.37 


1.45 


51.56 ±2.82 ±0.45 


0.15 


70.19 ±3.26 ±0.37 


1.55 


43.51 ±2.56 ±0.39 


0.25 


72.03 ±3.30 ±0.41 


1.65 


41.62 ±2.54 ±0.45 


0.35 


72.18 ±3.30 ±0.45 


1.75 


41.24 ±2.53 ±0.46 


0.45 


69.90 ±3.20 ±0.50 


1.85 


35.15 ±2.35 ±0.40 


0.55 


72.86 ±3.25 ±0.59 


1.95 


29.72 ±2.23 ±0.36 


0.65 


73.68 ±3.27 ±0.64 


2.05 


23.80 ±1.98 ±0.33 


0.75 


64.64 ± 3.09 ± 0.58 


2.15 


18.04 ±1.68 ±0.29 


0.85 


66.59 ±3.13 ±0.63 


2.25 


18.47 ±1.70 ±0.33 


0.95 


64.59 ±3.10 ±0.67 


2.35 


9.35 ±1.17 ±0.19 


1.05 


64.04 ± 3.09 ± 0.70 


2.45 


8.20 ±1.17 ±0.19 


1.15 


63.11 ±3.10 ±0.67 


2.55 


3.24 ±0.80 ±0.08 


1.25 


56.26 ±2.92 ±0.55 


2.65 


2.15 ±0.85 ±0.06 


1.35 


51.57 ±2.82 ±0.45 


2.75 


2.14 ±1.02 ±0.06 
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